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Abstract 

The next to leading order (NLO) quarkonium-hadron cross section is calculated in perturbative 
QCD. The corresponding leading order (LO) result was performed by Peskin more than 20 years 
ago using the operator product expansion (OPE). In this work, the calculation is performed using 
the Bethe-Salpeter amplitude and the factorization formula. The soft divergence appearing in the 
intermediate stages of the calculations are shown to vanish after adding all possible crossed terms, 
while the collinear divergences are eliminated by mass factorization. Applying the result to the 
Upsilon system, one finds that there are large higher order correction near the threshold. The 
relevance of the present result to the charmonium case is also discussed. 
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I. INTRODUCTION 



Twenty five years ago, Peskin and Bhanot and Peskin [2j showed that in the heavy 
quark limit, the interaction between a hadron and a heavy quark bound state could be 
described in perturbative QCD, and calculated the dissociation cross section of the quarko- 
nium by a hadron to leading order (LO). According to their calculation, which was based 
on operator product expansion (OPE), the J/ip dissociation cross section by a hadron was 
found to be very small and in the order of /ibarn in the threshold region. Later the result 
was re-derived and the target mass correction calculated 0, 0, IE 0; Du t the qualitative re- 
sult remained the same. Initially, such small cross section strongly supported that the J/ip 
suppression seen in relativistic heavy ion collision at SPS was a consequence of QGP (quark 
gluon plasma) formation^. On the other hand, other approaches, such as the quark ex- 
change mo del jl El , the meson exchange model [HI EH 12 .Hp . "LH laL the QCD sum rule 
method 17, ll8Ul£|. and other non perturbative methods 20| |. predicted a much larger cross 



section of few mbarn. The discrepancy is particulary large near the threshold region 21 
Although the model calculations themselves have large uncertainties and model dependen- 
cies, it is generally believed that such discrepancy exists because the QCD LO calculation, 
especially near threshold, is valid only for a very large quark mass, larger than that of the 
bottom quark. However, no systematic analysis has been worked out in this context, as the 
formalism based on the OPE is quite complicated even in the LO. 

A few years ago, Oh, Kim, and Lee [fjj used Bethe-Salpeter amplitude and factorization 
formula to reproduce Peskin's result on the dissociation cross section. Because this method 
is relatively simple, it opened the possibility to calculate the higher order correction, which 
will be carried out in the present work. 

There are two types of initial and final states in the NLO calculation. One is $ + q —>■ 
Q + Q + q, and the other is & + g — > Q + Q + g, where $ is a quarkonium. In the course of the 
calculations, collinear divergence, infrared divergence, and soft-collinear divergence appear. 
Infrared divergence disappear after adding the one loop diagram with Q + Q final state, while 
the collinear divergence is eliminated by mass factorization. Dimensional regularization in 
the MS scheme is used throughout this work, including the parton distribution function. 
The counting scheme, introduced by Peskin, to systematically study heavy bound states is 
used and applied to NLO. In addition, the large N c limit are taken throughout this work. 

In Section 2, the Bethe-salpeter amplitude and the LO calculation are reviewed. In 
Section 3, the elementary cross section of $ + g — > Q + Q + q is calculated and mass 
factorization is introduced. In Section 4, the calculation for the § + g^Q + Q + g process 
is presented. The effective four point vertex is introduced in section 5. It is then shown how 
the infrared divergences disappear when the relevant crossed terms of Born and one loop 
correction are included. Mass factorizations in the process <3> + g — ► Q + Q + g are presented 
separately when the emitted gluons are hard and soft in section 4 and 5 respectively. In 
Section 6, the result is applied to the upsilon dissociation cross section. Limitations when 
applied to the charmonium case is also discussed briefly. Appendix A summarizes the 2 and 
3-body phase space. Appendix B contains the derivation and the list of angular integration 
used in the present work. Appendix C gives detailed calculations of $ + g — > Q + Q + g 
diagrams. Appendix D gives some comments about order counting. 
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FIG. 1: The Bethe-Salpeter equation for quarkonium 



II. LO (<£ + g -> Q + Q) CALCULATION 

Here, the LO result, re-derived by OKLp using the Bethe-Salpeter amplitude, is pre- 
sented again for completeness. 

The Bethe-Salpeter equation, represented diagrammatically in Fig. (JTJ), is written as 

— — 7 Q zA(tf + Pl + p 2 )T fl (K + Pl + p 2 , K) 

xtA(K) la V(K + p 2 ), (1) 

where Cp = (N 2 — 1)/N C , iA(K) is the quark propagator, and iV(K) the gluon propagator. 

In the heavy quark limit, the Kq contour integration is dominated by the residue at 
ko = —m — K 2 /2m + ie over that at ko = —qo — m — K 2 /2m + is. In this limit, V(K + p 2 ) 
becomes three dimensional, 

V{K + P2 )\ Ko= _ m _ R2/2m ~ 1 (2) 

\K + p 2 \ 2 

If further, ^(g, p) is defined as, 

/ ^: iA ^ 1 ) r ^^ 1 ' _P2 ^ A ( _?j2 ) = ( 3 ) 

where q = p\ + p 2 , p = {p\ — p 2 )/2, the Bethe-Salpeter equation becomes, 



dp 



2tt 



d 3 K - , -t _ 1 



Ml,P) = l 9 °F / l^ A (Pi) / TTT^^M^K + q/2)l aT7 ^—- A{-p 2 ) 



(2tt)3' Hl ^' a \K-p\ 2 \ 



g 2 C F [d 3 Kl + 1 ° 1-7° 1 m 

-7 M<i> K n<* — — , (4) 



e +i?/mj (2tt)3 2 '™ ' ,a 2 |£_^s 

in the heavy quark limit and in the q rest frame. e G = 2m — m$ is the binding energy of 
the quarkonium. If we assume ^(q, p) to have a structure of ^2~lig l J—2~'4 ) {\p\) ■, the Bethe- 
Salpeter equation reduces to the nonrelativistic Schrodinger equation for the Coulombic 
bound state, 

(e + f/mmp\) = g 2 C F J ——^\K + p\), (5) 



whose spatial form is, 

2/i 47r r 
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if)(r) = ei/j(r), (6) 



where \i = m/2 is the reduced mass. 

At the same time, the Bethe-Salpeter vertex reduces to 



q q \ ( p 2 \ , , , /mi 1 + 7° * 1 — 7 



2 " 2 V \" m/™^ 17 V iV, 



r M = k + y^w^^- (7) 



* c 



Fig. (J2J is the leading order diagrams for quarkonium dissocation $ + g — > Q + Q- Among 
them, the first three diagrams are of the same order in the large N c limit, while the last two 
diagrams are suppressed by l/N c . Such N c countings are easily obtained using the double 
quark line representation for the gluon lines, since one notes that compared to the first 
two diagrams, the third diagram has additionally a quark loop and two factors of g. This 
altogether gives a factor of N c g 2 , which scales as 0(1) in the large N c limit, as the coupling 
g scales as l/(iV c )a. In contrast, the last two diagrams do not have an additional quark loop 
while they have the two additional factors of g, and hence are suppressed by 1/N C compared 
to the third diagram. Technically, this can be rephrased in terms of color matrices, as the 
last two diagrams carry the color factor T b T a T b = —T a j (2N C ) and hence is suppressed by 
1/N C compared to third diagram p. 

The counting scheme in scales starts from noting that the binding energy e a = 
m(N c g 2 /167r) 2 ~ 0(mgr 4 ). This suggests from Eq. (jSJ) that the three momenta of the heavy 
quarks are of 0(mg 2 ). The next important step is to take the external gluon momentum 
\k\ and its energy k° to be of 0(mg 4 ), which are smaller than the typical heavy quark mo- 
mentum of 0(mg 2 ) inside the bound state. This is the essence of the factorization in the 
present approach; namely, the separation scale is taken to be of 0(mg A ) so that the bound 
state property of 0(mg 2 ) can be taken into account as Wilson coefficients. Then, from the 
energy conservation m$ + k = 2m + \p 1 \ 2 /2m + \p" 2 \ 2 /2m, one has 

\ft\ ~ \p 2 \ ~ 0(mg 2 ), k° = |fc| ~ 0(mg 4 ). (8) 

Counting in the NLO process <& + q(g) — * Q + Q + q(g) are obtained similarly by assuming 
that the incoming and outgoing parton momentum are of 0(mg A ). That is, 

\fi\ ~ \P2\ ~ 0(mg 2 ), k° = \h\ ~ k° 2 = \k 2 \ ~ 0(mg A ), (9) 

where k\ is the incoming quark (gluon) momentum, and k 2 is the outgoing quark (gluon) 
momentum. Under this order counting scheme, quark propagators are expanded as below. 

1 + 7° 1 1 + 7° P-k * 
A(p + k) = — - — — H — — + 



2 k + is 2 m(k + is) 2 2mk + is 

1 + 7° -1 1+7° -P-k $ 
A{p — k) = — - — — H 7: Ti ~W + 



2 k — is 2 m(k — is) 2 -2mfc + is 



1-7° -1 I-7 -p-k 
A(-p + k) = —r- 1 —. - + 



2 kg — is 2 m(ko — is) 2 2mko — is 

A(- -k) = 1 ^'^ - § (10) 

{ P ' ' 2 k Q + ie 2 m{k + ie) 2 2mk + ie 1 ' 



4 




FIG. 2: Leading order diagrams for + g — > Q + Q. The two lower diagrams are suppressed in 
large N c limit. 



Here, p is the on-shell momentum of heavy quark. The first term is of l/mg 4 order, and the 
next two terms are 1/mg 2 order. The third diagram of Fig. (j2J seems to be higher order 
than previous two diagrams with respect to the coupling g. But they are of same order 
under above counting scheme. Detail is given in Appendix D. 

Using the Bethe-Salpeter amplitude Eq. (JJJ) and the heavy quark propagators Eq. (fTUjl. 
the leading order amplitude may be derived as, 



xn( Pl )^^7,^T^T^(p 2 ). (11) 



The total cross section then becomes 

3 

.7.2 ( A _ I 



where A = q ■ fci/m$, and the quarkonium is assumed to be in the coulombic IS state, 

VV^(p1 = ia|32Vi (M fl f +1)8 , (13) 

and a Q = 16tt / '(g 2 N c m) is the Bohr radius. The coupling g and the heavy quark mass m 
are determined by fitting the measured quarkonium spectrum, such as ip and ip' for the 
charmonium states, with those of the coulomb bound states. Details of the derivations are 
given in ref. pj. 

The hadronic cross section is obtained by folding the partonic cross section with the 
parton distribution function, using the factorization formula, 



0$ ^(A) 



f 1 

I dx(T^—p ar t on {xX 1 Q s )Dp ar t on {x ) Q\ (f4) 
Jo 
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FIG. 3: Next to leading order diagram with initial quark & + q — > Q + Q + q 



where x is the momentum fraction of a parton, and Q is the separation scale. As has been 
stated earlier, the scale Q is set to be the binding energy, which is of 0(mg 4 ), then it is 
natural to include the 'bound state' properties obtained from momentum scale of 0{mg 2 ) 
in the Wilson coefficient. 

At this point, it should be noted that the factorization theorem in Eq.()14|) is valid only 
when the mass of the quark is very large, so that the binding energy of 0(mg 4 ) becomes 
larger than the typical hard scale > 1 GeV. Otherwise, the two standard sets of corrections 
will not be small. Namely, the higher order correction to the perturbative cross section will 
not converge and the higher twist effects will not be negligible. 

If the heavy quark is not sufficiently heavy, the time scale involved in forming the bound 
state will not be short enough compared to the typical time scale involved for the bound 
state to interact with the external partons, and hence, the contributions from multiple 
gluonic interactions will not be negligible. Such multiple gluonic effects correspond to the 
higher twist effects. Moreover, even if one calculates the Wilson coefficients for such higher 
twist effects, nothing much is known about the higher twist distribution functions inside 
the hadrons, and the corresponding contribution to the hadronic cross sections can not be 
calculated. 

The perturbative cross section of the leading twist also has its own problems when the 
quark mass is not heavy enough. To begin with, to implement Eq. (jl4|) at the separation 
scale Q = e Q , one needs the parton distribution function D parton (x, e) defined at that scale. 
Moreover, the perturbative calculation for the leading order cross section cr$_ parton (xA, e) 
may not be convergent. Nonetheless, all such questions can be answered by the explicit 
NLO calculation, which will provide an quantitative estimate of the correction to the LO 
cross section, and also determine the valid energy range of the LO cross section. 

The binding energies for both the charmonium and the Upsilon systems are around 0.75 
GeV, and the aforementioned corrections are potentially not small. However, by applying 
our formal NLO calculation to the Upsilon system, we will investigate the convergence 
and the valid energy range of the leading twist cross section. As for the parton distribution 
function, we will use the MRST parton distribution function at its minimum scale Q 2 = 1.25 
GeV 2 , and investigate its uncertainty due to the variation in the scale to larger values. 
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III. NLO PROCESS <!> + q^Q + Q + q 

A. Collinear divergent elementary cross section 

Fig ((5J) represents the lowest order diagrams involving the quarks, § + q^Q + Q + q. 
The invariant matrix element for this process is given by 



M" NLO - q 



gu{k 2 )7 v T a u{k 1 ) 



x -9 



{h - k 2 ) ■ 



dip{p) 
dk 



g + (fcio - k 20 )^— g. 



-/ n 1 + 7o 



xu(p 



-li9 



fi% 



1 -7o 



2 2 
The averaged square in D-dimension is 



T a v(p 2 ). 



d Pj 



(15) 



7T71-2 2 12 4 2 / e D \ a fcio - k 20 - e c 

= T 779 m m *U (k 1Q - k 20 f 

2 2fci • k 2 1 J 4Jfei ■ k 2 



k 20 ) 2 ^ 



It may be modified to the following covariant form, 



(16) 



|M| 



NLO-q 



2 n 4 2 6 /e \2S + il-2m $ e c 

— 7T# m 2m $ — — 

3 vm/ (s + m) d 



x 



s 2 + w 2 
(2m$) 2 t 



p-4) 



yg +M) 

2t 



(17) 



where s = 2g • ki, u = —2q ■ k 2 , and t = —2k\ ■ k 2 . 

The parameterization of the three body phase space follows ref. (12]. The initial and final 
momenta are set to the following, 





= (E g ,0,. 


..,0,0, \p\ simp, \p\ cosy? - k w ) 


h 


= {kw,0, 


..,0,0,0, few) 


Pi 


= (Ei,0,. 


.., 0, —k 2Q sin 6\ sin 9 2 , —k 20 sin 9\ cos 9 2 , —k 2 o cos#i) 




= (E 2 ,0,. 


..,0,0, \p\ sinyj, \p\ cos </?) 


A- 2 


= (k 20 ,0, 


. . , , k 2 o sin #i sin 9 2 , /C20 sin 9icos9 2 , k 2 o cos 6*i ) . 



The following new invariant variables are introduced, 



ti = (ki - p 2 ) 2 - m 2 = -2ki ■ p 2 

u i = (q - P2) 2 -m\-m 2 = -2q ■ p 2 

s 4 = (k 2 + pi) 2 - m 2 = 2k 2 ■ pi = s + ti + Mi, (19) 

where s = (q + ki) 2 . 
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Eg, kio, Ei... may be expressed in terms of the invariant variables s(s), ti,ui, S4. 



E q 



E 2 
k 20 

\p\ 



s + ui + m% 



2^s A + 
s + ti — m% 
2 V Si + m 2 
s 4 + 2m 2 



2\/s4 + m 2 
ti + Mi + 2m 2 



2^j r s A ~^rri? 
s 4 



2\/s4 + m 2 

y/(u! + h) 2 - Am 2 s 



2 v / s 4 ~T 



txs 4 - s(ui + 2m 2 ) + (s 4 - s + 2m 2 )m% . . 

cosy? = . — . (20) 

(s + ti - m%) y/(ti + Ui) 2 - Am 2 s 

Using these relations, Eq. (|T7|) are expressed by five variables s(s),ti,ui,9 1 ,9 2 , because 



u = (q — k 2 ) 2 — m\ = 2k 20 (—E q + (\p\ cosy? — fcio) cos Q\ + \p\ sin ip sin 61 cos 6 2 ) 

{ = (fa - k 2 f = -2k 10 k 20 (l - cosfli). (21) 

The differential cross section for the three body decay is 

2 d 2 a _ 1 1 /i~ D+4 (su x t x -m\t 2 -m 2 s 2 \ V sf" 3 



dhdm 2(4n) D T(D-3)\ s 2 y? ) ( S4 + m 2)#-i 

x f de 1 sm D - 3 e 1 [ de 2 sm D - 4 e 2 \M\ 2 . (22) 

Jo Jo 

The derivation is given in Appendix A. When 9\ = 0, t = and the term 1/t in Eq. (|17|) 
gives collinear divergence. Defining as below, 

I&fi= d9 1 sm D - 3 9 1 d6 2 sm D - 4 6 2 -, , (23) 

Jo Jo 2 p{u + sy' 1 ; 

and expanding it with respect to D — 4 

I™> = ^-/'f + / (1J) + 0(D - 4), (24) 
the differential cross section of § + q^Q + Q + q is regularized as below, 



S 




2 faNLO-q _2 S 1 4 2 m N4fM f S 4 

_L_ ( _ /(^) + 2(s + m*e )I ( lf - 2s(s + 2m $ e )/^ 5) + 4s 2 m $ e /^ 6) ) 

-2m%I^) + 4m|e /(°' 6 ) - J^' 3 ) + 2{s + m*e )4 M) 
-2s(s + 2m*e )4 1,5) + 4s 2 m $ e 4 1,6) 

+ l 7 - B + ln ^ — 

X ( - J^ 3) + 2(s + m*e ) J^ 4 * - 2s(s + 2m$e )/^ 5) + 4s 2 m$e / a i 6) ) 
-^5 3) +m,e /5 4 ) ]. (25) 
From the list of integration in Appendix B, 

/(^) = *L I = 2n 1 (26) 

- 1 a (A + By a (sX)^ K ] 

where X = — [ui + m%)/(s + ti), which is 1 — k 2 o/k w in the quarkonium rest frame. The 
definition of a, A, and B are given in Appendix B. The terms proportional to l/(D — 4) 
come from collinear divergence, and are eliminated by mass factorization. 



B. Mass factorization 

Collinear divergence is eliminated by mass factorization, which moves the divergent con- 
tribution to the parton distribution function. When one parton is seen with finer scale, it 
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is not seen as one parton, but a sum of several collinear partons as shown in Fig. (|%|). In 
other words, parton distribution function has scale dependence, and the collinear parton 
with transverse momentum less than the scale should be included in the parton distribution 
function. Only the parton with larger transverse momentum is included in the perturba- 
tive calculations. Therefore collinear partons should be subtracted from the perturbative 
calculation. In the MS scheme, mass factorization is defined as 



, 2 d <j NLO - 

S : : 

dtidui 



, 2 d 2 (JNLO-i 



dt\du\ 2n 



a s / 1 dx 

— I — Pji(x) 







X 



+ 7s + In ■ 



D-A 



47T/T 



f d2&LO -\ (27) 



dt\du 



where s = xs, t\ = xt\. The parenthesis means the integration of transverse momentum from 
zero to momentum scale Q. aNLO-i is the next leading order cross section of quarkonium 
and parton i. &NLO-i is the reduced cross section after mass factorization, which is finite. 
o~LO~j is the leading order cross section of quarkonium and parton j, and crio-j = &LO-j- 
Pji is the splitting function of parton i to parton j. In & + q^Q + Q + q process, 



2N r 



x 



x 



(28) 



is needed. The LO differential cross section is 



g 2 C ^ L ,° 3 = ^-g 2 m 2 (2m^) 4 f-) 2 SX ?™* e ° 5(sx + t x x + u x + m%) 
dtidui 3A/ C \mJ (sxp 

D-2 1 / suiti - m%t\ - m 2 s 2 \ ^ 

X 2 T(D/2 - 1) V Anfi 2 s 2 ) ' 1 ' 

The LO differential cross section can be obtained from substituting Eq. (fTTj) into Eq. (jAlOj) 
in Appendix A. 

After mass factorization, the reduced differential cross section is 



, 2 d 2 CTNLO-q 

s — 

dt\du\ 



1 



fl V(2m$) 4 - 



e„\2 1 



3 (4tt) 2 ^ "" v \mJ S + t! 

- mKs + h) s 4 /(0;5) + 2m$eo/( o,6) : 
7r S4 + m 2 

I 1 ( i n s l + /(La)- 

{sXf V Q2( S4 + m 2) ^ o 



2(s + m<s,e c 
(sX) 4 



In 



6 4 



Q 2 (s 4: + m T 



(1,4)*\ J±_ 

+ I ° ) + (sxy 



2s(s + 2m$6„, . , 

+ (sxy ' ln 

Amq,e s 2 - , 
m 



Q 2 (s4 + m T 



+ P 



(1,5)* 



where I (1J) * = a(A + Byi^/ir = a(sXy J< 1J) /tt. 



(1,6)* 



- 1 



(30) 
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FIG. 5: Dalitz plot for $ + q -> Q + Q + g. 



One should also note that the cross sections to NLO in Eq. (|3()|) and to LO in Eq. (j!2|) 
have the same large N c scaling, as the coupling constant g scales as \/yW c . 

The threshold of differential cross section in (J3U|) is s = 2m$e„. This comes from the 
quarkonium wavefunction (j!3|) . However, the physical threshold is s = 2m$e + e 2 . The 
term e 2 was ignored because g is of 0{mg A ). We circumvent this problem by substituting 
e G + e 2 /2mq> for e G in the differential cross section. 

C. Dalitz plot 

Fig. (j3J) is the Dalitz plot, which is drawn under the following two conditions. 



S4 
|2 



cos x| 



s + ui + ti + m% > 
-2uis + (s + m%)(ti + Ui) 



(^1 + M i) 2 — 4sm 2 



< 1, 



(31) 
(32) 



here x is the angle between q and P2, or q and pi + ki in CM frame Eq. ()A8|) . The elementary 
total cross section is obtained by numerically integrating Eq. (|30[) over the Dalitz plot. 
Furthermore, the hadronic cross section is obtained by folding it with the quark distribution 
function. 



IV. HARD PART OF THE $ + g ^Q + Q + g PROCESS 

A. Soft and/or collinear divergent elementary cross section 

Fig. (jnj represent the diagrams for the process $+g — > Q + Q+g. Among them, diagrams 
(13), (14), and (15) are ignored because they are higher order in g compared to the rest of 
the diagrams, in the present counting scheme, where the momentum of the internal gluon, 
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< 1 > 




< 2 > 



< 3 > 





< 10 > 




< 4 > 



< 12 > pi 



> 

k2 q 



k2 
-p2 



klS 



< 5 > 




< 7 > 





FIG. 6: Next to leading order diagram with initial gluon § + g — > Q + Q + g. 

which binds heavy quark and antiquark, is of 0(mg 2 ) (from Eq. and that of the external 
gluon is of 0(mg 4 ) (from Eq. 0). Details are given in Appendix D. 
The invariant amplitude in quarkonium rest frame is 



A *-Hu\(a,b) _ 
m NLO-g - - 



The detailed derivation of this result is given in Appendix C. It was checked that the 
following current conservation conditions are satisfied. 

%M^ Q _ g = k lv M^ - 9 = hxM^ g = 0. 
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The averaged squared amplitude is 



-2 2 
\ M \nlo- 9 = —ng 4 m 2 (2m$ 



11 



,6 I 2 

rn 



1 



2m$e„ 



s + u 



(2m$) 2 t 



D-2/2S 2u u 



2 V m 
D - 2 / (s 2 + u 



'is 2u u 2 s 2 \ 

—+—+—+— 

u s s z u z / 



2m|t 



V su /J 



(34) 



The first line is of order t, the second line of order t°, and the third line of order t 1 . 
Differential cross section from the first line is 



3 (4tt) : 



e Q . 5 s 4 



where is defined as follows, 



m s 4 + m z 

/( -l,2) /( -l,3) /( -l,4) J( -1,5) X 2m # 6 ( _ 1;6) 



~J2~ ) 



(35) 



/*7T /*7T 

Jo 7o 

= d6 1 c 
Jo Jo 



d6 2 sin D - 4 6 2 



1 



(*)*(«)' 
sin D - 3 0! sin D - 4 fl 2 



a — a cos 6*i)*(A + £> cos #i + C 1 sin Q\ cos #2) j 



and, A = —2E q k 2 o, B — B, and C — C. The products of invariant variables are decomposed 
as below. 



(s + u) 5 u us 5 (s + -u)s 5 (s + -u) 2 s 4 (s + u) 3 s 3 (s + w) 4 s 2 



(s + m) 5 s 



(36) 



The first line has no divergent term. 
Differential cross section from the second line is 



4 2/ 

-g m 



(2m,f ( *L) 1 [ j(o,D ( 3i-8m«e N 



2 8 1 
3"(4^> 

D -2 / s Ay /suiti - m\t\ - rn^P x D ' 4 (02) t -s + 2m$e \ 
+ 2r(D-3)V 47rs/iVs 4 + m 2 / ' V s 4 / 



+/ (o,i) ^^3£+_8m$£o j | j - (0i 2)^ -2g + 6m $ e o ^ 
+/(0j3) ^-2£-F4m*e^ 



2i + 4m^ + 4m^, /(0i4) _ 2/(0>6) + 4m$eoJ (o,6) 



(37) 
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j(o,2) j s S oft divergent term, because it is proportional to 1/u 2 , which is 1 / ^20 m quarko- 
nium rest frame. 

Differential cross section from the third line is 

^ 4 m 2 (2m*) 4 (^ V ' 



3 (47r) 3 \m/ S4 + m 



,2 



x 



1 , SasJ su\ti — m'it 2 — m 2 s 2 1 
+ IE + In V , L - ^ =L= + - 



D-A ,1J An^ 2 s^s A + m 2 2 

27r f s — 2m$e„ 1 s — 2m$e D 1 2m$e„ 1 
x— -r, tt t, ^ + 



a V s 3 s(x — 1) s 3 sx s 2 (sx) 2 

-2 / + (s + 4m$e ) / - s(s + 2m$e ) \ + 2m$e s 2 — 1 



v sxj 3 (sx) 4 (sx) 5 (sx) 6 

£— 2m$£o (i,i) _ £— 2m$£o (1.1) 2m$e (i ;2 ) _ 9r (i,3) 
+ s '3 J o g , 3 + p J o 2i o 

+ (s + 4m $ e )4 1 ' 4) - s(s + 2m q> t )I { ^ ) + 2m <I ,e s 2 4 1 ' 6) ] , (38) 
where J^' is expanded with respect to D — 4 



= + + 0(D - 4), (39) 

and, 

m = 2 JL \ = 2 JL \ (40) 

a (A + By a (s(X-l)y 

1/(D — 4) term of Eq. (|38|) is collinear divergent. Additionally l/a(l — X) term gives 
soft divergence, because both a and 1 — X are proportional to &2o- Thus this term gives 
soft-collinear divergence. 

B. Mass factorization in hard gluon emitted region 

The Dalitz plot for $ + g — > Q + Q + g is the same as $ + g — > Q + Q + q. But it is 
separated as soft gluon emitting region and hard gluon emitting region as shown in Fig. (J7|). 
The boundary line is S4 = A which is an arbitrarily small value. That is, if S4 is smaller 
(larger) than A, it corresponds to the region where soft (hard) gluon are emitted. Because all 
infrared divergences exit in soft gluon emitting region, hard gluon region has only collinear 
divergence. This collinear divergence is eliminated by mass factorization. 

, 2 d 2 (JNLO-g 2^ 2,7 NLO-g 

s = s - 

dt\dui dt\dui 

a s f 1 dx / 2 Q 2 \ -2 d 2 a LO 
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where P gg {x) is a gluon to gluon splitting function, and may be separated into hard part 
and soft part, which are proportional to 6 function and 5 function respectively 

2 



P gg (x) = N c 6(l -x-5) (- — - + - - 4 + 2x - 2x 2 ) 
+ W (l-x)(21n ff + H 



11 iV^ 

= 6 (l - x - 5)P g H g + 8(1 - x)P g s g . (42) 

We ignore the factor proportional to iVy(the number of light quark flavor) in the soft 
part, because it is suppressed in the large N c limit. The boundary of hard and soft mass 
factorization is x = 1 — 5. 5 is related to A by 5 = A/(s + ti), because x < 1 — 5 means 
S4 > 5(s + ti). After hard part mass factorization, Eq. (}3*Hj) becomes 



e G \ 2 1 



3 (An) 2 y ' \mJ s+t x 



x 



s - 2m $ e Q 1 / sj (1,1)* 

P (l-X)s\ n Q 2 (s 4 + m 2 ) 
s-2m*e 1 / sj (1)1> 

s 3 sA V g 2 (s 4 + m 2 ) 
2m $ e 1 / s 2 (1>2) 



i n z± i_ n 

s 2 (sX) 2 \ Q 2 (s 4 + m 2 ) 

^_ ( i n ?1 + /(i- 3 )* 

sA) 3 V Q2( S4 + m 2) ^ o 



1 / s 2 
-(i + Am,e ) — (in + tf^ 

1 / s 2 
+ 2m ^o) ( In g2(g / +m2) + 4 1>5> 
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where J^* = a(A + B)J^' V /tt = a(X - 1) J^/tt. After mass factorization, collinear 
divergence l/(-D — 4) of Eq. (JHHj) is removed. 



— 2m$e„s 



r(i.*) 



1 



In 



sX) 6 V Q 2 (s 4 + m 2 



+ 4 



(1,6)* 



(43) 



V. SOFT PART OF THE § + g^Q + Q + g PROCESS 
A. Differential cross section for soft gluon emitted part 



In hard gluon emitted region, differential cross section is the sum of Eq. fl33j) . Eq. (|3*T|). 
and Eq. (j4*3^) . But in soft gluon emitted region, t — > 0, u — > 0, and s 4 — > 0. In this limit, only 
soft and soft-collinear divergent terms contribute and the differential cross section becomes 



d Z <7 



dt\du\ 



where 



Ofs + Ei+Mx) / as 4 s — — — 
y dt\du\ 

3 (47t)^ \m 



5(s + ii + mi) 



X 



1 A-i/sttiti — m|t? — m 2 s 2 1 
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D-A ' A7Tfi 2 ms 2 
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A a/ SMiti — m|r 2 — m 2 s 2 



+ ln z 



A a/ switi - m%t\ - m 2 s 2 
+2 In — - : — • In 



ms 



+(2 7 £ + 



47rm/x 2 s (s + ti)ra$ 

A a/ SMiti — m|r? — m 2 s 2 
1) In — — : 



+1 E + IE 



47T/i 2 m$(s + Ei, 
7T 2 

Y + 2 



(44) 



V A-VB 2 + C 2/ V ^4 + 5 



i _ y/B 2 + c 2 l 2 A + V& 2 + & 

+ ftr ; hr (45) 

A + B 2 A-VB 2 + C 2 

and we used the limiting values A 2 — B 2 — C 2 — > s 4 m|/m 2 , a(v4 + 5) — > ss 4 /(2m 2 ), and 
(A + _B) 2 /(A 2 — B 2 — C 2 ) — > [sm/m<f,(s + ti)] 2 . For the definition of soft gluon differential 
cross section, refer to j23|. In Eq. (jH|) , the first line in the square bracket comes from soft 
divergence, and the others come from soft-collinear divergence. These soft divergences may 
be eliminated by adding the mixed term of Born diagram and its one loop corrections. 
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B. Effective four point vertex 

Before considering the one loop correction, it is helpful to introduce the effective four 
point vertex. This vertex is attached to a quarkonium, a gluon, a heavy quark, and a heavy 
antiquark line and defined as 



\p\ 2 \ gVKg) vj 



x 



1+ 7o i » l-To T a 
2 r% 2 



(46) 



It is just the leading order invariant amplitude given in Eq. (jlljl . except that fcio is 
substituted by \p\ 2 /m + e . Although \p\ 2 /m + €q = k w to LO, it is not so in general. Using 
this effective vertex, the matrix element for the process & + g^Q + Q-\-g represented in 
Fig. (0) can be reproduced by five diagrams shown in Fig. (JHJ). Specifically, 



fif\(a,b) 



-gu( Pl ) 1 x T b A(p 1 +k 2 )M^ a >(k 1 ) 



(a), 



g 



+M^(k 1 )A(-p 2 - £; 2 ) 7 A T b \v(p 2 ) 

AO 



k 20 + is 



u(p 1 )M'»(k 1 )[T a ,T«\v(p 2 



(47) 



where M^ a >(&i) = M^ v (k l )T a . 

M 2 is the same as Mi with (/c 1; a, v) and (— k 2 , b, A) exchanged. 
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FIG. 9: One loop corrections. 



M 



,v0 



-k\o + is 



u( Pl )M^(-k 2 )[T b ,T a }v(p 2 ). 



M3 is a diagram which emits a gluon from the external gluon leg. 



M ^X(a,b) = _ ig fabc u ^ M ^(c)^ ki _ k ^ v ^ 2 y 



1 



(fa - A/2) 2 + i£ 
x Pi + k 2 ) a g vX + (fcx - 2k 2 y 9 x a + (-2k x + k 2 ) x gl] 



The sum of all diagrams is exactly the same as Eq. ([33)1 . 



M\ 



fiuX(a,b) 



NLO-g 



(48) 



(49) 



(50) 



Introduction of this effective vertex has some benefits. It makes the calculation much 
easier and one does not need to consider the inner structure of the four point vertex, which 
is very complicated when considering the one loop corrections. 

C. One loop correction 

Several comments are in order before considering the one loop correction. In dimensional 
regularization, 



T(D/2) f d D q A 



7T 



D/2, 



,D;-2a 



A 



D v -2a 



-q) a Dj/2 - a D v /2 



a 



(51) 
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Dj is the dimension which regularizes the infrared divergence, and Djj the ultraviolet diver- 
gence. A is the cutoff of the momentum integral. Generally this integration is zero. But we 
have left over a = 2 case, because it shows clearly the cancellation of infrared divergence and 
ultraviolet divergence separately. Second, in contrast to the order of typical loop momentum 
appearing in Eq. (jSJ), which is of 0(mg 2 ), we set the order of gluon loop momentum in the 
one loop corrections to be of 0{mg A ). This is to explicitly separate the soft part that cancels 
the soft divergence coming from emitted gluons of 0(mg i ). Keeping these comments, all 
divergences may be eliminated systematically. 

Diagram < 1 > of Fig. Q is the one loop correction of heavy quark and antiquark 
external lines. The product of on-shell heavy quark propagator and its self energy is 



iA( Pl )£(pi) = -ig 2 J (0j7 M A( Pl + kh,A( Pl )T a T a 1 



k 2 + ie 



9 n „l + 7o f d D k Am 2 1 1 . , 

y 2 J (2n) D ( Pl + k) 2 -m 2 + iep 2 -m 2 + iek 2 + ie' K ' 

We assumed that p\ is slightly off-shell, and used the heavy quark propagators as in 
Eq (JTUJ). Moreover, A(pi + k) is replaced by [A(pi + k) + A(pi — k)}/2. Then one has, 



iAtpOEfrx) = ig 2 T a T 



a l+ 7 o f d D k 



1 



(2tt) d (k + ie)(k - ie)(k 2 + ie) 



(47r)Vr( 



A 



D-A 



A 



D-A 



LD/-4 D 



u 



4J 



a 1 +70 



(53) 



In the ko contour integration, the residue at ko = ie or at ko = —ie makes the one loop 
correction pure imaginary. Therefore, in these cases, the mixed term of LO and its one loop 
correction vanishes. The self energy of the antiquark is the same except that (1 + 7q)/2 is 
replaced by (1 — 7o)/2. However it has the same contribution to the differential cross section, 
because spinor wavefunction u{p\) is proportional to (1 + 7 )/2, while v{p2) is to (1 — 7o)/2 
in the heavy quark limit. 

As can be seen from Eq. (J53)) . the renormalization constant of the heavy quark mass 
has no divergence. Only the renormalization constants of the heavy quark and antiquark 
fields have both ultraviolet and infrared divergences. The differential cross section from the 
mixed term of Born diagram and the same diagram but with the self energy insertion to the 
external heavy quark or antiquark line is 



2„BV1 



,d 2 a 



dtidui 



3 (An) 2 
x8(s + h + mi) 



4 2/o \Al € o 

g m (2m$) I — 



m 



s — 2m$€„ 



Dj-A D 



u 



Al 



(54) 



Note that the result was divided by 2, because half of the divergence is used for the 
renormalization constant of external quark antiquark wavefunction. It may be checked that 
the above infrared divergence cancels the soft divergence of Eq. (j44j) . 

Diagram < 2 > of Fig. Q is the external gluon line correction. The self energy of the 
external gluon is 
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ba 



d D k 



(2k 2 + 2k-k ± + 5k 2 )g au + (AD - 8)k a k„ 
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+ (2D - A)(k a k lv ) + (D — 6)ki a ki„ 
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:9 z N c 8 ba g, 



(k + k x ) 2 k 2 
(6 - 8/D)A; 2 



1 /• d D fe 
'"" ' 7 (27r)^L(P-x(l-a ; )(-^))2 



+ 



(2a; 2 - 2x + 5)£; 2 
(F-x(l-x)(-A; 2 )) 2 



(55) 



Here, we ignored quark loop contribution, as it is suppressed in the large N c limit. Terms 
proportional to k\ a and k\ v will vanish due to the current conservation condition of the LO 
amplitude. Assuming k\ to be slightly off-shell, the gluon self energy may be expanded with 
respect to k\. 



1 | 2x(l-x)(-k{) | 



(k 2 - x(l - x)(-kf)) 2 k 2 ' k 4 

1 1 2x(l-x)(-k 2 



(k 2 - x(l - x^-k*)) 2 
Keeping only the l/k A terms, 
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(57) 



The differential cross section obtained from multiplying this and the Born cross section 



is 



2^BV2 



,d 2 (T 



dt\du\ 3 (47r 



-g 4 m 2 (2m$ 



xS(s + t 1 +u 1 )- 





1 



§ s — 2m$6„ 



s' 4 



4 



(58) 



Again, the result was divided by 2 for the same reason. 

Diagram < 3 > of Fig. (JHJ) is the quarkonium external line correction. For quarkonium 
external line, there is no direct one loop correction in QCD. But we assume that its one 
loop correction is the same as that of heavy quark and antiquark lines in Eq. (J54)) . This 
assumption can be proven to be true from noting that the quarkonium is the bound state 
of a heavy quark and an antiquark, and therefore its field operator is a composite operator 
composed of a quark and an antiquark field. 

Diagram < 4 > of Fig. © is another type of one loop correction. However it vanishes 
and has no contribution. 
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(59) 



Finally diagram < 5 > of Fig. (JUJ) is the vertex correction. Soft collinear divergence of 
Eq. is eliminated by these diagrams. Such cancellation can be anticipated, because the 
soft collinear divergence of Eq. (J44)) and the cross section coming from multiplying diagram 
< 5 > with a Born diagram are coming from the different cutting of a two loop diagram for 
the self energy of the quarkonium with no infrared divergence. The amplitude becomes, 



M (a)V5 



2 rabc 



g 2 f 



d D k 

W) 



75«(Pi) W${k + h)A(-p 2 - k) lp T b 



x 



+ 1(3 T b A(p 1 + k)M$(k + k 1 )]v(p 2 ] 
-2h - kfgZ + (h - kY£ + (h + 2k) v g^ 



(k 2 + ie)((k 1 + k) 2 

1 i U 

"7 9i 



■ IE) 



g z^^ c u{ Pl ) l -±^^g» l —^T a v{p 2 ) 



d°- l k 



X 



2k l0 ^ ■ k l9u0 + k 10 (2k 10 



\k\)Wgl 



(27T) 

- 2k, 



D-l 



9ip(p) 



(k + hH 



Alk^iholkl + k-h) 
-2k l0 9 -^ ■ k l9vQ - k 10 (k 10 - \k\)^gi + 2(k - ■ kgi 



4|fc|(fci + \k\)(k 10 \k\ + k-ki 



(60) 



Two terms of the last equation are residues at ko = — \k\+ie, and at ko = —k w — \k+ki\+ie 
respectively. In the second term, k + k% was substituted by k. If each momentum is set to 
be 
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(61) 
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then, 
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(62) 



where means transverse direction with respect to k\. It is manifest that current con- 
servation condition k"Mffi V = is satisfied from Eq. (}62j) . because Eq. (}53j) and Eq. (|H7jl 
are intrinsically zero. The differential cross section from M y5 , after multiplying it with the 
Born diagram, is 
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(63) 



The double pole 1/(D — A) 2 is cancelled with that of Eq. 



D. Coupling constant renormalization and soft part mass factorization 

The ultraviolet divergence in one loop correction may be removed by renormalization of 
the coupling constant g. 



9b -> 9 



1 + ^ 



7i \D -A 
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Atx^ 



00 



(64) 



where the renormalization scale is set to the heavy quark mass. (3q — ^N c in the large N c 
limit. 

The sum of soft differential cross section Eq. (jUj) an d the terms obtained by multiplying 
the Born term and its one loop correction Eq. (|54j) . Eq. (J58)) . Eq. (|63|) still has collinear 
divergence, where Eq. (J54j) should be doubled because of quarkonium external leg correction. 
This remaining divergence is removed by the soft mass factorization in Eq. (J42)) . which 
corresponds to substituting the second part of Eq. (pEZJ) proportional to 5(1 — x) into Eq. ()41|). 
becomes divergence- free. 
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dtidui dtidui ° dtidui 

Below is the summary of the elementary total cross section for & + g^Q + Q + g. 
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Here, g 2 is the sum of Eq. (J35J), Eq. (|3*7jl . Eq. (jl3^1 . The first line and the second line 
depends on A (or 5). But their sum is independent of it, because it appears as the lower cut 
in the first line and as the upper cut in the second line. The results are shown in Fig. (fT0|) . 
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VI. UPSILON DISSOCIATION CROSS SECTION 

As an explicit application, the above result is applied to the upsilon dissociation cross 
section. The two independent parameters of the theory is determined by fitting the physical 
masses of mx(is) and m r (2S) to the energies of the Coulomb bound states. Specifically, from 
the relation m T (2S) — ^T(is) — 3/4e Q , the binding energy is found to be 750MeV. Also, 
the bottom quark mass is found to be 5.1 GeV from equating it to (mx(is) + e o)/2. The 
coupling constant g is then found to be 2.53 from g 2 = lQn /N c \Je /m. 

The left and right graphs in Fig. (fTU|) represent the elementary total cross sections of 
$ + g — > Q + Q + q and $ + g — > Q + Q + g respectively. In both graphs, there are 
regions of energy where the cross sections become negative. These negative cross sections 
originate from mass factorization, where finite parts of the differential cross section have 
been subtracted out and put into the definition of the distribution function. Therefore, 
the cross section becomes physical only after folding the elementary cross sections with the 
parton distribution functions (PDF) and adding them to the LO contribution. 

To obtain the total cross section, we used the MRST2001LO PDF for the LO re- 
sult, and the MRST2001NLO PDF @ for the NLO. We used the PDF calculated in the 
MS scheme, because our perturbative calculations, including the subtractions in the mass 
factorization, were performed in the MS scheme. If different schemes were used in the per- 
turbative calculation and in the PDF, the scheme dependent finite pieces would not match, 
and the result would be inconsistent. In the original scheme of Peskin, the scale of the PDF 
was be taken to be the binding energy of the system, which is 0.75 GeV in the present 
system. However, in the present example, we will take it to be 1.25 GeV 2 , which is the 
minimum Q 2 scale in the MRST PDF. 
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FIG. 10: The left figure is the elementary cross section for & + q^Q + Q + q and the right figure 
is that for $ + g — > Q + Q + [g). In the right figure, (a) the dashed line is the Born term given 
in Eq. (|12[). (b) The Dash dotted line is the hard gluon part, namely, the integration of the sum 
of Eq. Q35|) . Eq. (|37|). and Eq. (|43j) over the hard part of the Dalitz plot plus the integration of 
the 'In 5' dependent part of Eq. H65|) over the soft part of the Dalitz plot, (c) The Dotted line is 
the soft plus one loop correction, namely, the integration of Eq. (|65|) excluding the 'In <5' dependent 
part over the soft part of the Dalitz plot. The solid line is the sum of (a), (b), and (c). 

The left graph of Fig. (jllj) shows the total dissociation cross section to LO and to NLO. 
The ratio between the cross sections calculated to NLO and to LO, plotted in right graph 
of Fig. (jTTj) . shows that perturbative QCD approach is acceptable only in a limited region 
of energy and large corrections exist in the threshold region. 

The separation scale of 1.25 GeV 2 is low, making it questionable whether one can apply 
the present formalism to the Upsilon system. On the other hand, one can not take the 
separation scale to be arbitrarily large in this example, as one would invalidate all the 
counting schemes and the non relativistic approximations used in deriving the formula. 
Nevertheless, to asses the uncertainties associated with takin the scale of the PDF to be 
low, we modified the scale Q 2 of PDF to 2.0 GeV 2 , and compared the result to that obtained 
with Q 2 = 1.25 GeV 2 . As shown in Fig. (JED, the total cross section changes by less than 
10% for y/s < 25GeV . As shown in the right graph of Fig. (|12|). 25GeV is also the upper 
limit of the window of energy region where the ratio between NLO to LO is minimal. Hence, 
although the scale dependence is non negligible in the present example, the uncertainties 
are within the estimated errors coming from the perturbative expansion. 

We also applied the present calculation to the J/ip dissociation cross section. But in 
J /if) case, the soft plus one loop correction has large negative value making the hadronic 
cross section negative. This suggests that the formalism breaks down for the charmonium 
system. As the quarkonium is heavier, the relative contribution of this negative part is 
smaller. Therefore we conclude that charm quark is not heavy enough to use the present 
formalism in the present form. 
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FIG. 11: The left figure is the T(lS)+nucleon total cross section to LO (dashed line) and to NLO 
(solid line). The right figure is the corresponding ratio between NLO and LO results. 
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FIG. 12: The left figure is the T(lS)+nucleon total cross section with Q 2 = 1.25 GeV 2 (straight 
line) and with Q 2 = 2.0 GeV 2 (dashed line). The right figure is the ratio between NLO and LO 
when Q 2 = 2.0GeV 2 . 



VII. CONCLUSION 

We have reported on the NLO calculation for the quarkonium parton cross section in 
QCD. All the collinear divergences have been shown to cancel through mass factorization 
and the soft divergences among themselves. The result constitutes an exact QCD calculation 
at the NLO in the formal heavy quark limit. 

Explicit application to the Upsilon system shows that there are large NLO corrections 
especially near the threshold, as has been originally anticipated by PeskinQ. Neverthe- 
less, we have identified a window of energy range where the NLO are under control, such 
that the perturbative QCD results are reliable. Moreover, we have identified the origin of 
large corrections and assessed the uncertainties through the magnitude of the higher order 
correction. 

The application to the charmonium system confirmed that the discrepancies existing 
between LO QCD result and hadronic model result on the charmonium dissociation cross 
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section by hadrons especially near threshold are partly due to large higher order corrections 
in QCD. Nevertheless, since the separation scale is in the order of the binding energy, 
a thermal mass of few hundred MeV for the partons will be enough to soften the NLO 
correction and to make a perturbative treatment meaningful at finite temperature. 
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APPENDIX A 

Here the phase space of 2-body and 3-body decay are reviewed. The initial flux F is 



F = 4^{q ■ ki) 2 - m\m\ = 2s, (Al) 
and the total cross section of 2-body decay is 

/i- D + 4 f d D Pl f d D p 2 



2s 



[ d pi [ d p 2 x+ , 2 2\r+/ 2 2\ 

j j W) 5 ^ 5 {Pl ~ m)5 {P2 ' m) 



x{2tc) d 5 D { Pi + p 2 - q - h)\M\ 2 



2s 



I (|^ 5+ (P2-^ + ((? + ^i-P2) 2 -m 2 )|Mp. (A2) 



Because 



J d D p 2 5 + {pl -m 2 ) = ±J dp2od\p-2\ 2 \p-2\ D - 3 dn D - 2 5 + (pl - m 2 - m 2 ) 

J dE 2 dcosx(El -m 2 )^(l - cos 2 x)^ , (A3) 



D/2-1 



T(D/2- 1) 
where we have used the formula 



and 



dVL D _ 2 = sin D 3 #isin D 4 6 2 ... sin D _ z d9 1 d6 2 ---d9D-2 

2tt d / 2 - 1 

= nm^T) smD " 3hdh ' m 
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the cross section becomes 



o- = , /o ^n-^r/n/o TT / dE 2 d cos X (Ei - m 2 ) 2 (i_ cos 2 x ) 2 



2s (2n) D - 2 T(D/2 - 1) 



x5(s + t 1 + Ml )|M| 2 , (A6) 



where (? + fci — P2) 2 — m2 — s + t± + u±. 
In the center of mass frame of q and hi, 



q = (E q ,0,...,0,k w ) 
h = (k 10 ,0,...,0,-k 10 ) 

p 2 = (f? 2 ,0,...,|p|sinx,|p|cosx), (A7) 



and 



E q 



ho 
E 2 



s + m% 
2^ 
s 



2^ 
ti + Ui 

+ Mi) 2 - Asm 2 



COSY " 2MlS ~( S + + /AON 

sy (ti + U\Y — 4sm z 
From the above relations, the Jacobian from Ui,ti to _E 2 ,cosx is 

fs 

dE 2 d cos x = =dtidu 1 . (A9) 

Sy/(ti + Mi) 2 — 4sm 2 

Then the differential cross section is 



d 2 (j ii / 1 \ D / 2 rs-u^i — m|t 2 — m 2 s 2n ~ 



s 2 

dtidui r(D/2 - 1) V4tt7 



fJ, 2 S 2 



x5(s + t 1 + u 1 )\M\ 2 . (A10) 

Next, the 3-body phase space is 



/j ' 2(g " 4) /" ^2 f d D Pl f d D k 2 2 



x5 + (pl - m 2 )5 + {kl){2'K) D 5 D {q + h - Pl - p 2 - k 2 )\M\ 2 
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d D p 2 d D p5 + (p 2 2 - m 2 )5 D (q + k 1 -p 2 -p) 



1 /i~ 2 ( J - 4 ) 
~ 2i(2vr) 2 ^ 3 

x y d D k 2 d D p 1 5 + (p 2 1 -m 2 )5 + (k 2 2 )5 D (p-k 2 -p 1 )\M\ 2 . (All) 

In the above equation, new variable p is introduced and D- dimensional p integration and 
D-dimension delta function of p are inserted into phase space. The last line is the same as 
the two body phase space, and it becomes 

2_1HJ ™2\D-3 /-7T /-7T 



r(f - 1) - 



III 



2 r(D-3) (p2)#-i 



d8 1 sm D - 3 6 1 I sin 1 ^ 4 9 2 \M\\ (A12) 



o 



Moving to q and k\ center of mass frame, the differential cross section for three body 
decay is 



s 2 



d 2 a 1 1 /i~ D + 4 fsu l t l -m%t 2 -m 2 s 2 \ ! ^ ± sf" 3 



dhdm 2(An) D T(D-3)\ Pfi 2 J ( S4 + m 2)#-i 

x f de 1 sm D - 3 6 1 [ d9 2 sm D - 4 9 2 ~\M\ 2 . (A13) 
Jo Jo 

APPENDIX B 

Here the angular integration which are used throughout this work is derived. The required 
angular integration has the form 

= r d9 1 an D - 3 9 1 f d9 2 sm D ~ 4 9 2 -, 

' Jo ' tf(u+sy 

M r M sin^sin^ 

1 ' 2 <a + bcos9 1 ) i {A + B £0*9^0^9^0*9^' 1 ; 







Jo 



where 



A = 2E q (k 1Q - k 20 ) + 2k 10 (k 10 - \p\ cosy?) 

B = 2k 20 (\p\ cosy? - fcio) 

C = 2k 20 \p\ sin ip 

a = -2k w k 20 . 



(B2) 



D denote that it is D-dimensional angular integration. If b = —a, A 2 = B 2 + C 2 , the solution 
has an explicit form 2^1 23{. 



Aid) 



2ty 



F{D/2 - i - l)r(D/2 -j- 1)T(D - 3) 



(2aY(2Ay T 2 (D/2-l)T(D -i- j -2) 

r D A- B 
xF 1<2 1,3,-f- 1; 



2A 



(B3) 
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F12 is the hypergeometric function. 



F h2 (a, b, c; x) = — ^ / eft ^(l - t) M Q- - ^)~°, 

r(6)r(c- 6) i 



r(6)r( c 

which has the following properties 



Fi )2 (a, 6, c;x) = F lj2 (b, a, c; x) 
F 1>2 (0,b,c;x) = 1. 

In the present case b = —a, and A 2 ^ B 2 + C 2 . First, considering Ip' 1 ^, 



7T /*7T 

d#i / d0 2 
o Jo 

1 



sin D ~ 3 0i sin D - 4 2 



a(l - cos0i)(A + .Bcos^i + Csin0i cos0 2 ) 



a(A + 5)7 ^ (1-cos^) 



a(A + B)J Vo (1- cos 00 



1 



A + 5 



A + B cos 6»i + C sin X cos 2 



/1 + /2 



where / is separated into infinite part I\ and finite part I 2 in 4 dimension. 
From Eq. l|H3jt 



and 



r(l,0) 



7T 



A + 5 a(A + 4' 



1 



1 ^ J 

, d cos 0i / <i0 2 - — 

a^ + S)^ 7 (1- cos 0i 

A + B 



1 - 



A + B cos 0i + C sin 0! cos 2 
~ f 1 d cos 0i 



+ e(£>-4) 



a(A + 5) 1 - cos 0i 
A + B 



v /(A + 5cos0 1 ) 2 -C 2 sin 2 0i 



e(D-4) 



7T 



a(A + 5) 



In 



(A + 5) 2 - (1 - cos0i)(fi 2 + C 2 + AB) 

COS 0\ = + 1 



+ (A + 5) ^ (A + 5 cos 0i) 2 + C 2 sin 2 0! 

* m .^t B > 2 „ + e (0 -4)' 



COS 01=— 1 



e(D-4) 



a(A + 5) A 2 - B 2 - C 2 
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As a result, 



,-(1,1) 



7Y 



a(A + B) 



D — A 



+ In 



(A + Bf 



A 2 -B 2 -C 2 \ 



+ 0(D-4). 



(B9) 



Ip'^ with higher order j is derived by differentiating I^' 1 ^ with respect to A. 



" (j — 1)! dA?~^ ■ 



(BIO) 



Below are the summary of the results, (the subscript 'D' is omitted for simplicity). 



j(i,o) = 



2tt 1 
~~aD-4 



(Bll) 



/(i,i) 



7T 



a(A + 5) VD-4 
D-4 



+ In 



LA 2 - 5 2 - C 2 



+ 



+2Li 2 



A + VB 2 + C 2 



A + B 



2 A - VB 2 + C 2 

B+y/B 2 + C 2 \ (A-yfW+C* 

— -)- 2 H A+B 



A - VB 2 + C 2 
+0((D-4) 2 )) 



(B12) 



/(1,2) 



7T 



a(A + B)- 



(— 

\D- 



+ In 



(A + Bf 



A 2 -B 2 - C 2 



+ 



0{D-A)) 



(B13) 



j(l,3) = 



7T 



+ In 



a{A + Bf \D-A 

+(A + B) 2 



(A + B) 2 
A 2 -B 2 -C 2 \ 
A 2 + B 2 + C 2 
(A 2 -B 2 - c 2 y 



+ 



2A(A + B) 



A 2 -B 2 
3 + 0(D- 



-C 2 
4)) 



(B14) 



/(1,4) 



7T 



+ In 



(A + Bf 



+ 



2A{A + B) 



+{a + b) \a 2 

-11/3 + 0(D- 



A 2 + B 2 + C 



A 2 -B 2 -C 2 \ A 2 -B 2 -C 2 
2A 3 + 6AB 2 + 



-B 2 
4) 



C 2 )- 



+ (A + 5)^ 



3(A 2 - 5 2 - C 2 ) J 



(B15) 
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/(1.5) 



7T 



a(A + B) 5 \D - A 



+ In 



(A + Bf 



A 2 



A 2 + B 2 + C 2 



C 2 ) 2 
-C 2 ) 



+{A + B) 2 

+{ * A + B) \(A 2 -B 2 -C 2 Y 
25/4 + 0(^-4)' 



{A 2 - B 2 - 
AA 2 (B 2 



B 2 



C 2 \ 
B) 



2A(A + B) 

+ A 2_ B 2_ C 2 

3 2A(A 2 + 3B 2 + 3C 2 ) 



3(A 2 -B 2 - c 2 y 



1 



(A 2 -B 2 - C 2 f 



(B16) 



j(1.6) 



7T 



o(A + S) B VD-4 



+ In 



{A + Bf 



2A(A + B) 



A 2 _j_ r2 I /i2 



A 2 -5 2 -C 2 J A 2 -5 2 -C 2 

2A(A 2 + 35 2 + 3C 2 ) 



3(A 2 -B 2 - c 2 y 



+{A + Bf 
+(A + Bf 



{a 2 -b 2 - c 2 y 

AA 2 (B 2 + C 2 ) 1 
(A 2 -B 2 - C 2 Y + {A 2 -B 2 - C 2 ) 2 

32A 5 2A{3A 2 + B 2 + C 2 )^ 



(A 2 -B 2 - C 2 ) 4 



j(o,j) anc [ /( ij) ma y j-, e obtained from [22j and by the same method, 



5(A 2 -B 2 - C 2 ) 5 
137/30 + 0(D-A) 



(B17) 



j(o,U 



V^ 2 + c 2 



(B18) 



/(0,2) 



2tt 



A 2 - B 2 - C 2 



(B19) 



j(0,3) 



2ttA 



(A 2 - 5 2 - C 2 ) 2 



(B20) 



/(0,4) 



2tt(3A 2 + B 2 + C 2 
3{A 2 -B 2 - C 2 f 



(B21) 



j(0,5) 



2ttA(A 2 + E 2 + C 2 
(A 2 - B 2 - C 2 ) A 



(B22) 



7(0,6) 



2\2^ 



2tt(5A 2 + 10A 2 (5 2 + C 2 ) + (£ 2 + C 2 ) 
5(A 2 -B 2 - C 2 f 



(B23) 
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na 



A + ^/B 2 + C 2 



2B B 2 + C 2 + AB , 
+ -^t^ In 



(B24) 



j(-i,2) 



7ra 



2(AB + 5 2 + C 2 



(B 2 + C 2 )(A 2 - B 2 - C 2 ) 
B 



A + VW+C 2 

(B 2 + c 2 f/ 2 U - vs 2 + c' 2 



In 



(B25) 



j(-l,3) 



27ra(A + 5) 
(A 2 - E 2 - C 2 ) 2 



(B26) 



( _ 14) _ 27ra(3A 2 + 4A£ + B 2 + C 2 
~ 3(A 2 - B 2 - C 2 f 



(B27) 



/(-1,5) 



27ra(3A 3 + 5A 2 B + (3 A + B)(B 2 + C 2 )) 
3{A 2 -B 2 - C 2 ) A 



(B28) 



/(-1,6) 



27ra(5A 4 + WA 3 B + (10 A 2 + 6AB)(B 2 + C 2 ) + (-B 2 + C 2 ) 2 ) 
5{A 2 -B 2 - C 2 ) 5 



(B29) 



APPENDIX C 

Here the detail calculation of Fig. © is presented. Mi through M 8 correspond to tree 
diagrams (1) to (8). 
Mi is 



M 



(i) 



{ig) 2 u(p 1 )T^(p 1 ,p 1 - q)iA(pi - q)i„T a 



xtA{-p 2 - k 2 ) lx T b v{p 2 ) 
1x9 



2 m $- ( \ 



1+7° 1-7° , 



K 2 o 

n n P2 ■ fc 2 „o P2jffl nPlifli 

~9u9\ zr^2 y^zn h #a 



20 



m/c 2 o m ^ 



20 



(CI) 
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The first line in the bracket is of l/(mg 4 ) order, and the second line is of l/(mg 2 ) order 
from Eq. Q. As will be shown later, the sum of I /(rug 4 ) order terms from Mi to Mg 
vanishes. As a result, the l/(mg 2 ) order becomes the leading order. 

M 2 is 



Kux = (ig)*u(p 1 )j x T b iA(p 1 +k 2 ) lv T a iA(q-p 2 ) 



2 m ^- ( \ 



1 + 7 C 



-liQ 



1-7°, 



T b T a v(p 2 ) 



+9u9x 



1 X mk 2( 
oo-l%) 



20 



[-ki + h) 



(C2) 



1 ho dp 

The third line in the bracket comes from the expansion of the Bethe-Salpeter amplitude. 



r M (pi + k 2 , -p 2 + k x ) = r M (pl + k 2 , f + k 2 ) 

= r M (pi,pi) + — — k 2 



(C3) 



r M (-P2, <? - P2) = r M (pi - k x + k 2 , fx -kx + k 2 ) 

= ^(pi.pi) H ^ {-kx + k 2 ) 



(C4) 



r M (pi - fci, -P2 - fe) = r M (pi - - fci 



r M (pi,pi) + — ^ 



(C5) 



M 3 



is 



M, 



(3) 



(ig) 2 u(px)j\T b iA(px + k^T^px + k 2 , -p 2 + kx) 
xiA(-p 2 + kx)i v T a v(p 2 ) 



x 



a7 



0„0 



u(px) 



1 + 7° 



"7*# 



7° , 



9 V 9\ 



+ \f\ 2 /m 



k w k 



10^20 



+9°A 



f -kx P2-kx- fx ■ k 2 f -k 2 , 

+ — + — — )^(P) 



mk 2 



mkink 



10^20 



mk 2 



+ [9i9x 
+9°A 



oP2j kxo — k 2Q g jPxj kxo — k 20 \ 



m k 10 k 20 

( 2 pi • fc 2 



m fci fc 2 o 



mkink 



10^20 



^10^20 dp 2 
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1 + 7° 



9 \ljf u {Pi)^—^- ( l 



1 11 - 



7 



T b T a v{p 2 ) 



x 



Q k w — k 2Q 



10^20 



P~2 • &i Pi-h- Pi ■ h Pi ■ k 2 



mk 2 



mk w k 20 



m k 2 o 



9 v 9x 



m k w k 20 



. o o ( fi • h + pi ■ fc 2 , , 



^20 . o ^10 — k 20 

+ 9u9\ 



m k w k 20 

kio ~ k 20 dipjp) 
k w k 20 dp 



■h 



(C6) 



Notice that |p2| 2 « \p\\ 2 — 2p~[ ■ (k\ — k 2 ) from the order counting of Eq. (jHJ). The last line 
in the last equation comes from the expansion of the Bethe-Salpeter amplitude in Eq. (|C3J) . 
Similarly, M 4 to M 8 are given as below. 



M, 



(4) 



{ig) 2 u(p 1 )^ u T a i/S.(p 1 - kx)Y^{pi - fci, -p 2 - k 2 ) 



xiA{-p 2 - k 2 ) lx T b v{p 2 ) 



2 /m # _, a+7 u 
5 \/^«(Pi)— 2~ 7*P 



7 



T a T b v(p 2 ) 



x 



■ o o~ e + |pil 2 / m ,/ \ 

9v9\ — t_ TO 



10^20 



Pi ■ h fi ■ k\ - p 2 ■ k 2 p 2 -k 2 . 

+ — — V(P) 



mk\ 



+ 9l9\ 



+9°A 
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mk 10 k 2Q 
P2j k 



m k 2 Q 



m k 



Kon 111 Ki 



> 20 



10^20 
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2pl • ki , . fcio - &20 <9V>(p) 



m$ _ . . 1 + 



7 



10^20 
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-li9 



-T b T a v(p 2 ) 



20 



o o ~~ k 

9u9x T T 

^10^20 



+9°A 



+ (glg x 
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mkf 



mkink 



10^20 



mk 20 



oPij k\o — k 20 q jP 2 jki — k 20 , 



m /cio^o 



m fci fc 2 o 



o o /Pi ' ^2 + Pi ■ h k 10 - k 20 dip(p) r 
-9u9x -t—t„ i>{p) + — ■ h 



mk w k 



20 



10^20 
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M, 



(5) 



(ig) 2 u(p 1 )T fl (p 1 ,p 1 - q)iA(p 1 - q)<y x T b 
xi^(-p 2 + k l ) lv T a v{p 2 ) 



■2 /-- ( .^l + 7°„,^l^/ :r6rMp2) 



g~\I^T u (Pi) 



9u9\ 



10 



I no g2vfcl nPljffA , oP2j^ 
v X ~mk 2 ~ h ^ A 
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10 



10 



(ig) 2 u( Pl )^T a iA( Pl - fc 1 ) 7A T 6 iA( ? -p 2 ) 
xr M (g-p 2 ,-P2)w(p 2 ) 



1 + 7° 



2 /"I*-/ x 

gj^ufa) 2 



-Tiff 



1-7° fc 



X 



1 



ff.ffAT-V'b) 
K 10 



J„0 Plj 



0jP2j „0„oPl " ^1 
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10 



'A; 
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^£a = igu(piW,(p 1 ,p 1 -q)tA(p 1 -q) 1 d T c v(p 2 ) 
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Xgfabciih + fc 2 ) 
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2/c 2 + A;i)^A5 + (-2fci + k 2 ) x g u ^j 



2 x 



i + tV^i-t 



-Tiff 
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x 
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ip(p) 



M S\ = igu(p!)^ d T c iA(q - p 2 )T^(q - p 2 , -p 2 )v(p 2 )- 



,,Srpc 



(h - k 2 f 
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1 + 7° 7 



-Tiff 
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Vi ■ h + p 2 ■ k 2 2kijP2k + 2k lk p 2j j A , , , 

9v\ + 9v9\ Hip) 

m m J 

+ {{ki + k 2 ) g uX + (-2k 2 + h) u g X o + (-2*4 + k 2 ) x g u oj 
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x ■ 
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Next considering others, 



(ig) 



d 4 k 
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-T b T a v(p 2 ) 
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xgf cdb ({2k + 2p 2 + k 2 ) x g aS + (fe - k-p 2 ) a g xs + (-A; -p 2 - 2k 2 ) s g a \) 

2 f™*.-.,„ A 1 +7° , ia 1 -'frrar 



-ng llx: 
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d A k 
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Because M^ 13) ~ are higher order in the order counting of Eq. (JHJ) compared to 

the other terms, they were ignored. Summing all amplitudes from AfW to M^ 16 \ the total 
invariant amplitude is given as, 



V op 



K 20 K l " K 2 



MUo^oV" 
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fcio fci • fa 



2 xl + 70 



10 



-7i5 



A" 



00 0J 

1 - 7o 



&20 M ' k 2 

[T\T b ] v(p 2 ). 



(C17) 



37 



APPENDIX D 



In this appendix, we summarize the present counting scheme and give details on deter- 
mining the scaling properties of certain diagrams. To begin with, the binding energy of the 
quarkonium scales as, 

e = m(N c g 2 /lQn) 2 ~ 0(mg 4 ). (Dl) 

In the quarkonium rest frame, the energy conservation condition for the process $ + 
q, 9(h) -> Q(pi) + Q(p 2 ) + q, g(h) is, 



i , \pi\ 2 \P2\ 2 mo x 

where eo is the binding energy of the quarkonium, and k w , fc 20 are the energies of the 
incoming and outgoing partons respectively. p[ and pi are respectively the three momenta 
of the heavy quark and antiquark from the quarkonium. From this relation, the following 
order counting can be deduced. 

\fi\ ~ \P2\ ~ 0(mg 2 ), k w ~ fc 20 ~ 0(mg 4 ). (D3) 

The counting for the internal gluon loop momentum K, which connects the heavy quark 
and antiquark within the bound state, can be deduced from Eq. (J3j). Since the left hand 
side of Eq. (jSJ is of 0(mg 4 ), K must be of 0(mg 2 ). In contrast, the order of gluon momenta 
appearing in the perturbative one loop correction should be of 0{mg 4 ). This is so because 
the separation scale, which sets the cut off in the perturbative diagrams, are set to the 
binding energy, which is of 0(mg 4 ). Within the bound state loop, the internal energy and 
heavy quark propagator can be counted as 0{mg 4 ). 

From the above considerations, the order of each Feynman rules can be deduced and 
the results are listed in table (pj. Bound gluon means that it is the internal gluon, which 
produces the coulomb bound state, and whose momentum K scales as 0(mg 2 ). There are 
two types of three gluon vertex. In the first one, the vertex combines two bound gluons and 
one external gluon, while in the other, it combines three external gluons. 

The order of a diagram can be deduced from the above order counting scheme. For 
example, the left and the right diagrams of the Bethe-Salpeter equation in Fig.(Q) can be 
shown to be of the same order. Compared to the left diagram, the right diagram has an 
addition internal loop ((mg 4 ) x (mg 2 ) 3 ), two heavy quark propagator ((mg 4 )~ 2 ) , a bound 
gluon propagator ((mg 2 )" 2 ), and two coupling constant (g 2 ), which altogether gives order 
1, as the left diagram. 

The suppression of diagrams (13), (14) and (15) to the other diagrams in Fig. © may 
also be explained. Comparing diagrams (13) or (14) to (9), diagram (9) has additionally 
a heavy quark propagator {(mg 4 ) -1 ) and a quark gluon vertex (g), while diagram (13) or 
(14) has additionally a bound gluon propagator {(mg 2 )~ 2 ) and a three gluon vertex (two 
bound gluon + one external gluon of (0(mg 3 )). Hence diagram (13) or (14) is suppressed 
by 0(g 2 ) relative to diagram (9). Similarly, diagram (15) does not have a three gluon vertex 
(0(mg 3 )) nor a heavy quark propagator ((mg 4 )" 1 ) but an additional four gluon vertex (g 2 ) 
compared to diagram (9), and hence is relatively suppressed by a factor of g 2 . 

In certain cases, the simple counting scheme has to be implemented with care. As an 
example, consider comparing the order of the first two diagrams to the third diagram in 
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Feynman diagram order reference 

heavy quark (antiquark) propag ator (mg 4 )- 1 Eq. ifTHjl 

bound gluon propagator (mg 2 )~ 2 Eq. |(SJ) 

external gluon momentum nig 4 Eq. © 

bound gluon momentum mg 2 
three gluon vertex (two bound and one external gluons) mg 3 

three gluon vertex (three external gluons) mg 5 



TABLE I: Order counting for various Feynman diagrams 



Fig. (J2J). Using our naive counting scheme above, the third diagram can be shown to be 
suppressed by g 2 compared to the first two diagrams. However, in this case, the first two 
diagrams cancel to leading order in the counting and combine to give the same order as the 
third diagram. 
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